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ABSTRACT 

Aims: We present the detailed analysis of triple system KR Com with different observational tech- 
niques - photometry, interferometry, and period variation. 

Methods: The use of BVR photometry of the close-contact binary KR Com, which is the primary 
component of a triple system, helps us to better describe the properties of the components. The 
interferometric data obtained during the last 30 years sufficiently determine the visual orbit, but the 
use of minima timings of KR Com for the study of period variation together with the visual orbit is 
a novel approach in this system. 

Results: Basic physical parameters resulting from the light curve analysis agree well with the previous 
results from spectroscopy. The temperatures for the primary and secondary component resulted in 
5549 and 6072 K, respectively, and the amount of the third light in all filters is about 1/3 of the total 
luminosity. The distant third component revolves around the common barycenter on 11 yr orbit with 
a very high eccentricity (0.934) and this movement is also detectable via the period variation, which is 
clearly visible in the O — C diagram of times of minima observations. The use of minima times for the 
combined analysis helps us to independently determine the distance to the system (64.02 ± 9.42 pc) 
and also to confirm the orientation of the orbit in space. 

Conclusions: New minima observations and also spectroscopy would be very profitable, especially 
during the next periastron passage in the year 2017. 

Subject headings: binaries: eclipsing - binaries: visual - stars: fundamental parameters - stars: indi- 
vidual: KR Com 



1. INTRODUCTION 

Multiple stellar systems (i.e. of a multiplicity of three 
and higher) are excellent objects to be studied. Be- 
sides the statistics and their relative frequency among 
the stars, it is important to investigate these systems in 
detail, also because of their stellar evolution, their ori- 
gin, to test the influence of the distant components to 
the close pair, etc. 

One of these systems is KR Com (HD 115955, 
HIP 65069), which is also A component of the visual 
binary A 2166 (WDS J13202+1747). Its relative bright- 
ness is about 7.2 mag in V filter and its combined spec- 
tral type has been classified as F8V, lAbt I (| 19811) . The 
system has been discovered as a p hotome t rically variable 
one from the Hipparcos data by iSelam I (|2004[ ) . Its or- 
bital period is about 0.408 days only, but both eclipses 
are rather shallow (because of a close bright compan- 
ion). Although both minima are of similar depths, the 
system has been incorre ctly classified as a j3 Lyrae one, 
iKazarovets et al. I (|1999f ). 

The close binary s ystem has been extens ively studied 
spectroscopicallv bv lRucinski et all ()2002[ ). The radial 
velocity curve of KR Com has been derived and also an 
influence of the third body in the system to the close pair 
and its pro perties has been discuss ed in subsequent pa- 
pers by e.g. iD'Angelo et al. I (|2006l ). A very low value of 



the mass ratio of this system has been found, about 0.09 
only, which makes KR Com an exceptional case among 
the W UMa-type binaries. The photometry and light 
curve solution of the close pair is still missing, which led 
us to observe this interesting target. 

On the other hand, the close visual companion B (with 
an ang ular distance about onl y 0.1") was detected in 
1980 bv lMcAlister et al. I (|1983t) . The movement around 
a barycenter is very fast, its orbital period is about 11 yr 
and the eccentricity is very high about 0.9 . The orbit 
has been derived by lHartkopf et alH (|1996l ). Owing to 
the small angular separation between both components, 
any photometric and also spectroscopic observations can- 
not be done only for KR Com, and the influence of the 
third component has to be considered in the analysis. 
The eclipsing binary is the br i ghter component (Am~ 
0.6 mag ), and iRucinski et al. I (2002) give the luminos- 
ity ratio Ls/(Li + Lq) = 0.56. They also presented a 
broadening function, w hich indicated a s l owly rotating 
third component, and in ID'Angelo et all (|2006| ) the au- 
thors provided an estimate about the third-body mass 
M3 = 1.19 Mq, while the individual masses of the close 
binary are 1.42 and 0.129 M Q . 

Thanks to the unusually low mass ratio of the eclips- 
ing pair and the presence of the close interferometric 
component B we deal with a unique system, which is 
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Fig. 1. — Observations of KR Com from Hipparcos (plus) and 
ASAS (circles from 2009, dots from 2005, and triangles from 2003). 
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Fig. 2. — Long-term variation of brightness of KR Com (in 
quadratures). Dots represent the various catalog data, the plus 
signs stand for the ASAS data, and the open circles mark our 
observations. 

very useful to be studied in detail. 



2. PHOTOMETRY 

Owing to its relatively high brightness, the system has 
not been observed photometrically and studied in detail. 
There exist the light curves from the Hipparcos data (H p 
filter), and from the ASAS survey (V filter), see Fig. 
[TJ Because of the low photometric amplitude and high 
scatter, these can hardly be used for any analysis of the 
light curve. 

Thanks to the observations obtained by ASAS, which 
cover more than six years, we have found that there is a 
long-term photometric decrease of KR Com. It can be 
seen in Fig. [TJ where we present the light curves from 
from different epochs. As one can see, a steady decrease 
of its brightness is on the order of 0.1 mag during the 
last six years. 

The nature of this long-term behavior is questionable, 
as is which component of the system yields this varia- 
tion. On the other hand, we have tried to find some 
historical data in various catalogs, which contain at least 
some information about its magnitude in photometric fil- 
ters. VIZIER 1 lists 145 catalogs containing several data 
about KR Com, from which a few usable data points 
were selected. The most complete are the measurements 
in the V filter. In Fig. [2] we have plotted the data after 
1990. The long-term variation is clearly visible, neverthe- 
less this plot has a few flaws. For some of the points the 
photometric filter was not exactly the Johnson's V filer, 
but rather some "similar" one. And secondly, for some 
of the data points it was quite hard to find an exact date 
of the observation, only a mean epoch is presented. The 

1 http:/ /vizier. u-strasbg.fr/viz-bin/VizicR 



most deviating point from 1996 is a typical one. This par- 
ticular data point comes from "The PM2000 Bordeaux 
proper motion catalogue" (Ducourant et al. 2006), and 
in its description is mentioned that the filter used is a 
non-standard one, a combination of two different filters. 
Exactly the same problem arises for the second deviat- 
ing point from 1999, which comes from the "M2000: An 
astrometric catalog in the Bordeaux Carte du Ciel" (Ra- 
paport et al. 2001). However, both these magnitudes 
use the same combination of filters and the difference 
between them during three years is easily visible. A pe- 
riod of this variation could only be vaguely estimated, 
but it is probably longer than the period of the third 
body. 

One can speculate about the magnetic activity cycles 
in the system. Because of the spect ral type (F8V), 
these can play a role, see lHall I (|1989D . These cycles 
could slightly modulate the orbital period of the close 
pair. More o ver, a ccepting this hypothesis, according to 
lApplegatel (|1992f ) there is also a modulation of the lu- 
minosity of the star on the order of 5L/L ~ 0.1. If we 
propose that the most luminous component (the primary 
one) undergoes this variation, a variability of the bright- 
ness of the whole system results in only about 0.05 mag, 
which is much lower than observed since 1990. On the 
other hand, removing these two questionable data points 
discussed above, one can obtain a variation with an am- 
plitude of about only 0.1 mag, which could be described 
by this mechanism. Therefore, the nature of this varia- 
tion still remains an open question. 

We have observed the system during the seasons 2009 
and 2010. In total there are 17 nights of observations, 
but for the light-curve analysis we used only 7 nights 
of observations obtained from March 2010 to April 2010 
and carried out with the same telescope and detector at 
the private observatory by one of the authors (RU). Ow- 
ing to high brightness of the target, the refractor with 
a diameter of only 75 mm was used, equipped with the 
G2/KAF 0402ME CCD camera and standard B, V, and 
R filters according to specification by Bessell (1990). All 
the measurements were processed by the software C- 
Munipack 2 , which is based on aperture photometry and 
using the standard DaoPhot routines (Tody 1993). 

The observations were transformed into the standard 
system using the well-known transformation equations 3 . 
Star HD 115981 was used as a comparison, while the 
check star to control the non-variability of the two was 
star HD 116206. The atmospheric extinction has been 
neglected because the stars are very close to each other 
(11 arcmin), their spectral types are also similar, and 
the observations have never been obtained below 25 deg 
above horizon. We used this approach because our ob- 
servations have an accuracy not better than 0.01 mag 
(for the scatter in individual filters see below), therefore 
any additional correction on the order of 0.001 mag is 
practically useless for the transformation. The individ- 
ual exposure times range from 25 to 140 seconds. 

The rest of observations were used only to determine 
the minima times for a period analysis. These new 
ones as well as the already published ones are given 
in Table [TJ where the type of minima refers to the 

2 see http : //c — munipack.sourceforge.net/ 

3 see |http: / /brucegary.net / AUSky/x.htm| 
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TABLE 1 
Minima times of KR Com. 



HJD - 2400000 Error Type Filter Reference" 



48500.2095 


0.0075 


Pri 


Hp 


HIP 


48500.4200 


0.0129 


Sec 


Hp 


HIP 


52638.8769 


0.0165 


Sec 


V 


ASAS - [1] 


52758.41911 


0.00232 


Sec 


V 


ASAS 


52758.62161 


0.00301 


Pri 


V 


ASAS 


53058.4796 


0.0002 


Pri 


c 


Kraici (2005) 


53102.11786 


0.00255 


Pri 


V 


ASAS 


53150.2695 


0.0005 


Pri 


c 


Kraici ('2005') 


53475.0061 


0.0004 


Pri 


V 


Naeai, K. (20061 


53487.86332 


0.00518 


Sec 


V 


ASAS 


53488.05344 


0.00767 


Pri 


V 


ASAS 


53829.11608 


0.00468 


Sec 


V 


ASAS 


54225.87197 


0.00849 


Sec 


V 


ASAS 


54226.07038 


0.00326 


Sec 


V 


ASAS 


54233.4095 


0.0009 


Pri 


c 


[2] 


54580.40196 


0.00273 


Sec 


V 


ASAS 


54613.0406 


0.0002 


Sec 


Ic 


Naeai, K. (20091 


54921.47502 


0.00150 


Sec 


13 


this paper 


54921.47435 


0.00246 


Sec 


V 


this paper 


54921.47702 


0.00178 


Sec 


R 


this paper 


54921.47534 


0.00051 


Sec 


1 


this paper 


54924.53515 


0.00043 


Pri 


B 


this paper 


54924.53553 


0.00048 


Pri 


V 


this paper 


54924.53559 


0.00227 


Pri 


Pi 


this paper 


54924.53549 


0.00031 


Pri 


1 


this paper 


54940.44566 


0.00113 


Pri 


R 


this paper 


54971.44977 


0.00098 


Pri 


B 


this paper 


54971.45562 


0.00085 


Pri 


V 


this paper 


54971.45419 


0.00052 


Pri 


R 


this paper 


54944.52261 


0.00457 


Pri 


V 


ASAS 


54944.32228 


0.00565 


Sec 


V 


ASAS 


55259.67912 


0.00319 


Sec 


R 


this paper 


55274.37115 


0.00175 


Sec 


B 


this paper 


55274.36049 


0.00135 


Sec 


1 


this paper 


55274.36870 


0.00543 


Sec 


R 


this paper 


55274.35921 


0.00329 


Sec 


V 


this paper 


55274.56764 


0.00335 


Pri 


B 


this paper 


55274.56959 


0.00143 


Pri 


1 


this paper 


55274.56684 


0.00293 


Pri 


R 


this paper 


55274.56834 


0.00264 


Pri 


V 


this paper 


55280.49301 


0.00093 


Sec 


B 


this paper 


55280.49220 


0.00069 


Sec 


R 


this paper 


55280.48942 


0.00109 


Sec 


V 


this paper 


55281.50792 


0.00089 


Pri 


B 


this paper 


55281.50389 


0.00085 


Pri 


R 


this paper 


55281.50817 


0.00050 


Pri 


V 


this paper 


55293.53729 


0.00061 


Sec 


B 


this paper 


55293.53830 


0.00157 


Sec 


R 


this paper 


55293.53559 


0.00099 


Sec 


V 


this paper 


55294.56261 


0.00110 


Pri 


B 


this paper 


55294.56418 


0.00157 


Pri 


R 


this paper 


55294.56473 


0.00120 


Pri 


V 


this paper 


55357.38912 


0.00235 


Pri 


1 


this paper 


55358.41029 


0.00257 


Sec 


I 


this paper 



following ephemeris: 2454924.538 + 0.4079711 • E. For 
all these observations, the Kwee - van Woerden method 
(hereafter KW, Kwee & van Woerden 1956) was used 
for determining the time of minimum. This method is 
suitable for symmetric minima (which is the case for 
KR Com), because its principle is based on comparing 
the ascending and descending branch of the minima. 



TABLE 2 

Light curve parameters of KR Com. 



Parameter 


Value 


Parameter 


Value 


Tt [K] 


5549 ± 244 


Li (B) [%] 


50.67 ± 0.20 


T 2 [K] 


6072 ± 270 


L 2 (B) [%] 


17.75 ± 0.19 


q (= Ma/ATi) 


0.091 (fixed) 


L 3 (B) [%] 


31.58 ± 0.20 


e 


(fixed) 


Li (V) [%] 


47.63 ± 0.24 


i [deg] 


52.14 ± 0.46 


L 2 (V) [%] 


18.98 ± 0.23 


xi (B) 


0.774 


L 3 (V) [%] 


33.39 ± 0.24 


X2 (B) 


0.699 


Li (R) [%] 


47.05 ± 0.72 


x t (V) 


0.635 


L 2 (R) [%] 


19.04 ± 0.67 


X 2 (V) 


0.565 


L 3 (R) [%] 


33.90 ± 0.72 


x x (R) 


0.548 


Derived quantities: 


x 2 (R) 


0.485 


Ri [Rq] 


1.33 ± 0.04 


9l 


0.323 ± 0.014 


R 2 [R ] 


0.49 ± 0.02 


92 


0.351 ± 0.021 


JWbol.l [mag] 


4.34 ± 0.18 


Ai 


0.499 ± 0.027 


-Mbol,2 [mag] 


6.13 ± 0.24 


A 2 


0.484 ± 0.015 


a [mag] (B) 


0.00745 


Fx 


1.690 ± 0.066 


a [mag] (V) 


0.01583 


F 2 


1.647 ± 0.148 


a [mag] (R.) 


0.06737 



We do believe that our error estimates are much more 
reliable than those already published because first, 
we did a detailed analysis of the KW result for each 
minimum, which means that we used different data sets 
(neglecting some of the observed points) and compared 
these results. And secondly, we also used a polynomial 
fitting to determine the time of minimum and compared 
this result with the KW one. The error of the particular 
minimum light was computed as a maximum difference 
between all the different results from different methods 
from the mean value. All minima times given in Table 
[T] are heliocentric ones. Eleven new minima times 
were also derived from the data of the ASAS survey 
(Pojmanski 2002). 



3. THE LIGHT CURVE ANALYSIS 

For the light-curve analysis we used the program 
PHOEBE (Prsa & Zwitter 2005), which is based on the 
Wilson-Devinney algorithm (Wilson & Devinney 1971). 
The model atmospheres made by Kurucz (1996) were 
used. The derived quantities are: the individual tem- 
peratures Ti and T2, the inclination i, the luminosities 
Li, the gravity darkening coefficients gi, the limb dark- 
ening coefficients Xi, the albedo coefficients Ai, and the 
synchronicity parameters Fi. The limb darkening was 
approximated via linear cosine law, and the values of x\ 
and xi were i nterpo lated from the van Hamme's tables, 
Ivan Hammel (| 19931 ). 

Because the sys tem has been a n alyzed spectroscopi- 
cally in detail by iRucinski et al. I ([2002), we used the 
mass ratio value q — M2/M1 — 0.091 derived from spec- 




-0.6 -0.4 -0.2 0.2 0.4 0.6 

Phase 



Fig. 3. — Our new observations of the light curves of KR Com 
in three filters. The B and R curves are shifted for better clarity. 
The theoretical fits and parameters are described in the text and 
Table 



a HIP - Hipparcos observations, ASAS - data from the ASAS 716 
survey, [1] - see http:/ /var. astro.cz/ocgate l [2] - unpublished 
yet, see http://eclipsingbinary.web.fc2.com/vsnetmin.htm □> 12 
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troscopy for the analysis of our B, V, and R photom- 
etry. The other relevant parameters of the light curve 
were computed (with hxed circular orbit) and the results 
are presented in Table [2j The final fits and the photo- 
metric data are plotted in Fig. [3l The scatter in the B 
and V filters is about 0.010-0.015 mag, in the R about 
0.025 mag, and the numbers of data points are 451, 468, 
and 474, for B, V, and R, respectively. The individ- 
ual observational errors of the measurements were about 
5-10 times lower than a typical peak-to-peak scatter of 
the measurements in the particular filter. The same ap- 
plies for the measurement errors for the comparison and 
check stars. Their respective errors are also on the order 
of 0.003 mag only, but these values represent merely the 
mathematical results of errors from the photometry of 
the stars. More plausible physical errors are much above 
these values, about 0.01 mag. This is actually a very sim- 
ilar situation as for the minimum estimate and the error 
of the KW method as discussed in the previous section. 
No additional short-time variation has been found in the 
data. The BVR photometry used for the analysis is also 
available as an online table. 

The errors of the individual parameters were derived 
from the covariance matrix, following a standard proce- 
dure of error estimates, see e .g. "Numerical reci pes. The 
art of scientific computing" . iPress et al. I (|1986l ). 

From the light-curve analysis we see that the two com- 
ponents are in contact (A- type contact system), but their 
individual temperatures differ about more than 500 K. 
On the other hand, the difference between the tempera- 
tures of 523 ± 364 K is only marginally significant. We 
tried to find an appropriate solution with equal temper- 
atures (which resulted in a value of about 6100 K), but 
this solution resulted in \ 2 about 13% worse than the 
solution presented in Table [2J Another approach was to 
fix the temperature of the primary to T\ = 6100 K and 
to fit only the secondary temperature. This resulted in 
fit with x 2 about 9% worse than our final solution. Both 
alternatives have been tested via an F-test to gauge how 
liable they are to introduce another free parameter (tem- 
perature) to the model. We concluded that introducing 
this new parameter is significant at a level about 86% 
and 84.5% for these two different approaches. 

Our resulting values of parameters from the light curve 
fit can be compared with the estimates published in 
iRucinski et al. I (|2002f ). A value of the third light (33% 
of total light) agrees well with their estimate (about 56% 
of L\ + L2). The individual components are probably 
of F+G spectral types for the primary and secondary, 
respectively. The secondary component has apparently 
undergone a mass transfer to the primary. 

We were able to partly reveal the nature of the third 
body in the system via its temperature. This can be 
estimated thanks to the individual luminosities in the 
filters, from which one can derive the magnitudes in the 
particular filters. These values resulted in B = 8.96 mag, 

V = 8.39 mag, and R = 8.05 mag. As an independent 
comparison one can also check the value from the 
Washington Double Star Catalog (WDS) 4 , where the 

V magnitudes 7.78 and 8.38 mag are given for the pri- 
mary and secondary component of the visual pair, which 
excellently agrees with our value of 8.39 mag. From 

4 http://ad.usno.navy.mil/wds/ 



TABLE 3 

Interferometric observations of A-B pair. 



Date 


e 




p 


References 




fdeel 


[arcsec] 




1980.1566 


0.9 


U. 


1 A K 
. 140 


McAlister et al. (1983) 


1980.4817 


2.0 


n 
U. 


1 AO 
. 14z 


McAlister et al. (1983) 


1983.0510 


10.0 


n 
U 


117 
11 I 


McAlister et al. (1987) 


1983.0701 


12.2 


n 
u. 


1UO 


McAlister et al. (1987) 


1984.0558 


22.6 


n 
u. 


uoy 


HartkoDf et al. ( 1996) 


1986.4067 


327.7 


n 
U. 


UOO 


McAlister et al. (1989) 


1987.1194 


345.5 


n 
u. 


U f O 


Fu et al. (1997) 


1987.2642 


340.2 


n 
u. 


UOO 


McAlister et al. (1989) 


1987.2859 


343.0 


n 
u 


.uoy 


Fu et al. (1997) 


1987.2886 


348.8 


n 
u 


uoy 


Fu et al. (1997) 


1987.2914 


346.8 


n 


uoy 


Fu et al. (1997) 


1988.1022 


345.5 


n 
u. 


uy / 


Fu et al. (1997) 


1989.2300 


353.6 


n 
U 


11 ( 


McAlister et al. (1990") 


1990.2081 


357.9 


n 
U 


1 00 
loz 


Baleea et al. (1994) 


1990.2621 


356.0 


n 


-LOU 


HartkoDf et al. (1992) 


1990.2759 


356.0 


0. 


135 


HartkoDf et al. ( 1992) 


1991.25 


357.0 


0. 


136 


Perrvman & ESA (1997) 


1991.3187 


0.7 


0. 


140 


HartkoDf et al. ( 1994) 


1992.3098 


4.3 


0. 


138 


HartkoDf et al. ( 1994) 


1992.3127 


4.2 


0. 


139 


HartkoDf et al. ( 1994) 


1993.1971 


7.7 


0. 


128 


Hartkopf et al. (1994) 


2002.3224 


359.8 


0. 


140 


Horch et al. (2008) 


2002.3224 


0.0 


0. 


140 


Horch et al. (2008) 


2002.3224 


0.0 


0. 


139 


Horch et al. (2008) 


2006.1917 


17.8 


0. 


.077 


Mason et al. (2009) 



our estimates of the BVR magnitudes, the photometric 
indices B — V = 0.57 mag and V — R = 0.34 mag can 
be compared with the color- temperature relations, e.g. 
by iHoudashelt et al. I (|2000D . This gives temperature 
estimate of the third component of 5900 ± 200 K, which 
indicates a spectral type between F8 and G6. On the 
other hand, if we try to derive also the radius of this 
component from the relation between the radii, lumi- 
nosities and temperatures i?3 = ^/L^/Li ■ (T1/T3) 4 • R\ , 
we obtain values from 0.93 to 0.99 R©, which indicates 
a later spectral type, about G8 - Kl (Harmanec 1988). 



4. COMBINED SOLUTION FOR THE VISUAL ORBIT AND 
PERIOD VARIATIONS 

As we mentioned above, the A-B system revolves 
around a barycenter with a period about 11 yr. Because 
this orbit has been derived about 15 years ago and many 
new observations are available since then, there is a need 
for new up-to-date solution for this orbit. 

We used all available interferometric observations for 
the analysis of the pair collected in the "Fourth Cata- 
log of Interferometric Measurements of Binary Stars 5 " , 
lHartkopf et aT~l (|2001l ). There are 25 usable measure- 
ments (see Table [3]) . In some of these observations the 
position angle 9 has to be changed by 180° (i.e. quad- 
rant change, interchange of the components), the values 
presented in Table [3] are the corrected ones. 

A movement of the contact binary around the barycen- 
ter with the distant component has to produce a periodic 

5 http:/ /ad. usno.navy.mil/wds/int4. html 



The triple system KR Com 



5 





TABLE 4 


Final parameters of the long orbit. 


P 2ir £Liii6t 6r 


Value 


HJDo 


2453058 4fi40 + D 0021 


P Irlavl 
r |u<iyj 


n 40797003 + n 00000939 


P6 L ua J J 


4011.0 ± 98.4 


fvrl 


10.98 ± 0.27 


To 
1 u 


2442055.1 ± 118.0 


a [arcsec] 


0.1131 ± 0.0189 


cjr> [detrl 


301.8 ± 3.8 


e3 




i [deg] 


67.68 ± 2.03 


n [deg] 


210.03 ± 1.25 


A [day] 


0.0171 ± 0.0019 


/(Ma) [M ] 


0.3262 ± 0.0896 


M 3 [M Q ] 


1.598 ± 0.423 


oia [AU] 


3.677 ± 0.821 


a 3 [AU] 


3.563 ± 1.027 


7r [mas] 


15.62 ± 1.80 



2452000 2454000 
2000 2002 2004 2006 2008 2010 



0.14 - 



0.1 



0.04 - 




—0.02 



-0.06 -0.04 



-0.02 0.02 0.04 0.06 0.08 
RA[arcsec] 



Fig. 4. — Orbit of A-B pair on a plane of the sky. The dashed 
line represents the line of the apsides, while the dotted one stands 
for the line of the nodes. The cross indicates the position of the 
eclipsing binary. 

variation of orbital period of KR Com, a wel l- know n 
"light-time effect" (hereafter LITE), see llrwin I (|1959t ). 
The amplitude of the LITE depends on the inclination 
of the 11 yr orbit, the individual masses and also on 
the distance of KR Com. The distance from Hipparcos 
was originally derived as 76.51 ± 5.46 pc, i.e. parallax 
7T = 13.07 ± 0.87 mas (Perryman & ESA 1997), but 
the more recent value is 83.97 ± 5.3 2 pc, i .e. parallax 
7r = 11.91 ± 0.71 mas, Ivan Leeuwenl (|2007l ). Neverthe- 
less, the value of the parallax from the Hipparcos satel- 
lite could also be influenced by the movement on the 
long orbit. During the run of the satellite, the change 
of the position angle of the distant component was only 
about 15°, because in that time the B component was 
near apastron on its orbit. On the other hand, if one 
compares the error of its parallax with the other stars 
from the Hipparcos catalog with the similar parallaxes, 
the error is apparently not deviating significantly. 
For an analysis of the period changes of KR Com we 
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Fig. 5. — O — C diagram of KR Com. The primary minima are 
plotted as dots, the secondary ones as open circles. It is obvious 
that the respective errors of the published minima times should be 
larger than presented. 

collected all available minima observations. The very 
first ones are from Hipparcos data (recalculated primary 
minimum and derived also a secondary one), while the 
most recent ones are our new observations of 2010. All 
these heliocentric minima times are given in Tabled] We 
found that all the published minima times have their 
respective errors strongly underestimated (10 times or 
even more). 

We used a simil ar approach as described in 
IZasche fc Wolf I (|2007l ) for analyzing the interferomet- 
ric data together with the times of minima. The 
only difference was a calculation of both amplitudes - 
the semimajor axis of the visual orbit and the semi- 
amplitude of LITE in the O — C diagram. The main 
advantage of this approach is that we can indepen- 
dently calculate the distance (or parallax) to the sys- 
tem and compare it to the Hipparcos one. There- 
fore, the set of parameters to be computed is the fol- 
lowing: H JDQ,P,p3,TQ,i,a 1 uj 1 £l,e'$,A ) where a denotes 
the semimajor axis in arcseconds, while A stands for 
the semiamplitude of LITE. We used the least-squares 
method and the simplex algorithm (see e.g. Kallrath & 
Linnell 1987). One can also discuss the use of quadratic 
ephemeris because of some mass transfer between the 
components, but this was not used in our analysis due 
to insufficient coverage of the O — C diagram with data 
points in longer time scales. The two minima times 
based on Hipparcos data were only roughly estimated 
and therefore their use for the analysis is questionable. 

The analysis results in a set of parameters presented 
in Table IU where we list all computed values and some 
of the derived quantities. The final fits are presented in 
Figs. 2] and [5] For plotting the O — C diagram, the in- 
dividual minima in different filters were averaged to one 
point for better clarity. Figure [6] shows the residuals af- 
ter subtraction of the fit. As one can see, some long-term 
variation cannot be ruled out with the current data. This 
could be caused e.g. by a mass transfer between the com- 
ponents, but only future observations can confirm or rule 
out this possibility. Another potential explanation are 
magnetic cycles with a period of about 20 years (Apple- 
gate 1992), which cause not only a brightness variation, 
but also this period variation. A combination of both 
mass transfer and magnetic cycles is also possible. 

The distance to the system was derived from the fit, 
resulting in d = 64.02 ± 9.42 pc, closer than from the 
Hipparcos data. However, this value is only roughly es- 
timated because it is derived from comparing the values 
of A and a, but the detection of the LITE variation and 
deriving its amplitude is still problematic, and the 11 yr 
period is not yet sufficiently covered by data. Due to this 
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Fig. 6. — O — C diagram of residuals of KR, Com. 

reason, and due to the relatively large difference between 
our derived parallax and the value from Hipparcos, we in- 
troduced Fig. [7] In this figure we fixed the semimajor 
axis of the visual orbit (this value is much better defined 
than the amplitude of LITE) and plotted the values of 
the parallax tt and mass of the third component M3 with 
respect to the amplitude A of LITE. As one can see, our 
resulting value of A = 0.0171 days yields the values of 
M3 and 7r given in Table 21 If we consider the range of 
potential values of A according to its error, the values 
of 7r are still far from the value derived by Hipparcos. 
Moreover, the mass of the third body with a decreasing 
value of parallax is growing rapidly. On the other hand, 
if we presume that the mass of the third component is 
lower, about e.g. 1.1 M , then the amplitude of LITE 
is only about 0.013 days (which cannot reliably describe 
the fit in the O — C diagram) and the parallax results 
in approximately 16.5 mas, which is even more distant 
than the Hipparcos value. 

The other quantity, which can be compared with 
the previous resu lts, is e.g. the mass o f the distant 
component M3. iD'Angelo et al. I (|2006l ) presented a 
value about 1.19 M Q , while our estimation is a bit 
higher, about 1.60 M Q . This value has been derived 
from the assumption that the masses of the individual 
components are Mi = 1.42 M Q and M 2 = 0.129 M Q , 
which have been taken as constant input values. The 
resulting value of M3 is surprising, especially if we take 
into consideration that the luminosities in all filters 
resulted in slightly lower values (B: 46%, V: 50%, 
R: 51 % of Li + L 2 ) than assumed by ID'Angelo et al. I 
(2006). If the mass of the third body is higher and 
the luminosity lower, this could indicate that this 
component is also a binary, or the star is underluminous 
than one would expect for a star of the main sequence. 



5. DISCUSSION AND CONCLUSIONS 

The triple system KR Com revealed some of its basic 
physical properties for the first time. The light curve 
with shallow minima (due to a bright third component) 
shows that the system has a relatively low inclination of 
about only 52°. That we know the inclinations of both 
orbits is one aspect of KR Com system's uniqueness. 

On one hand, a difference between the inclination of 
the close pair and the inclination of the orbit of the third 
distant component (13 = 67.8°) is still relatively small, 
therefore one could speculate about a common origin of 
all three components (see e.g. Zakirov 2008). On the 
other hand, a presence of two orbits which are not strictly 
coplanar is a necessary condition for a so-called Kozai 
cycles (Kozai 1962), which could play a role in the for- 
mation of this system. As mentioned bv lD'Angelo et al~l 
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Fig. 7. — Diagram of semiamplitude of LITE versus mass of the 
third body and parallax of the system with the fixed semimajor 
axis of the visual orbit. The dash-dotted line represents the value 
of A from our final solution. 



(2006), the typical product of the Kozai mechanism is a 
close pair and a distant component on a non-coplanar 
orbit, which is the case here. 

The other computable quantity in this triple system is 
its nodal period, or the precession of the orbits, which 
can change the inclination of the close eclipsing pair and 
theref ore also the depths of the minima, see e.g. iMaver I 
(2005). Regrettably, the nodal period here results in a 
value of about 100000 yr. 

Today there are only 34 eclipsing binary systems 
known among the visual double stars that have known 
visual orbits. In all these systems one should expect a 
period variation of times of minima observed for these 
systems, but surprisingly, the LITE that agrees with the 
visual orbit has been detected in only seven systems (i 
Boo, VW Cep, C Phe, V819 Her, V772 Her, QS Aql, 
V505 Sgr). The system KR Com seems to be another 
example. The visual orbit with a period about 11 yr is 
well-covered by interferometric data, and the combined 
solution of visual orbit and O — C diagram of minima 
times is also dominated by the visual orbit. Therefore, 
the very first minima estimates based on Hipparcos data 
with their large errors were almost useless for the anal- 
ysis. The other interesting fact is that the already pub- 
lished minima show much larger residuals than their pub- 
lished errors. 

However, the combined fit brings new light into the 
system and confirms for example the u> and values. 
Because the same fit to the interferometric data could 
be realized by changing both angles by 180° (orientation 
of the orbit towards the observer), the LITE fit helps us 
to confirm the orientation of the orbit in the space with- 
out any need for spectroscopy. Nevertheless, the spectro- 
scopic data and new times of minima would be also very 
profitable, especially in the upcoming periastron passage, 
which will occur in 2017. 
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